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ABSTRACT 

Twin peak high frequency quasi-periodic oscillations (HF QPOs) observed in the power spec- 
tra of Low Mass X-ray Binaries (LMXBs), with either a black hole or a neutron star, have 
central frequencies that are typical of the orbital motion time-scale close to the compact ob- 
ject. Thus, twin HF QPOs might carry the fingerprint of physical effects in a strongly curved 
space-time. We study the azimuth phase <p{f) for orbital motion in the Schwarzschild metric 
and calculate the power spectra to check whether they display the features seen in the observed 
ones. We show that the timing of (pit) on non-closed orbits can account for the observed twin 
peak HF QPOs. The uppermost couple of peaks in frequency has the lower peak that corre- 
sponds to the azimuthal frequency v^, the upper one to + v r . 

The azimuth phase temporal behavior 0(f) on a slightly eccentric orbit in the Schwarzschild 
metric is described by a linear function of slope plus an oscillating term at the relativistic 
radial frequency v r . We deduce that the twin peak HF QPOs might originate from a frequency 
modulated (FM) signal driven by the kinematics of orbital motion in a curved space-time. 

Key words: X-rays: binaries - accretion, accretion disks - black hole physics - gravitation - 
methods: numerical 
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1 INTRODUCTION 

High freque ncy quasi-periodic oscill ations (HF QPOs) were dis- 
covered by Ivan der Kris et alj Jl99fjh . They are elusive modula- 
tions hidden i n the noisy X-ray flux from low mass X-ray bina- 
ries (LMXBs: iBelloni et alJl20l3) and are detected as Lorentzian- 
like peaks in their power spectra. Both neutron star (N S) and black 
hole (BH) LMXBs display HF QPOs (for reviews see Ivan der Klisl 
l2005tlRemillard & McClintockll2006h . HF QPOs often show up as 
couple, named twin peak HF QPOs. Because their time scale is 
typical of the orbital motion within 10 rJ3] from the compact ob- 
ject (milliseconds time-scale), twin peak HF QPOs might carry the 
fingerprints of the general relativity in a strong field regime. 

Twin HF QPOs are characterized by their central frequency, 
amplitude (root mean square; rms) and coherence (2-factor). In NS 
LMXBs extensive studies of both the rms and Q were pursued and 
revealed a s ystematic behavi or as function of the central frequency 
of the peak iBarret et al.l2006h . This was interpreted as possible due 
to the approa ch of the phenomenon to the innermost sta ble circular 
orbit (ISCO; iBarret et alj|2005t) . iKluzniak et ail H990t) described 
what type of fundamental information could be extracted from the 
fingerprint of the orbital motion on ISCO, such as constrains to the 
equation of state (EOS) of the matter at supra-nuclear density in 
systems with a neutron star. 



Twin HF QPOs are linked to the behavi or of the X-ray en- 
ergy spectrum of the source (Bellon3|20T3). It is known that 
the HF QPOs display correlations with the low -frequency QPOs 
dFord & van der Klisl 1 9981 ; iRemillard et alj|2002h . 

Despite their discovery since long, the puzzling twin peak HF 
QPOs still lack of an accepted modeling. A first m odel for the 
twin p eaks observed in NS LMXBs was proposed by Mill er et all 
dl998h . based on beat-fre quency mechanism s. In the relativistic 
precession model (RPM; I Stella et al] 1 19991) the twin peaks are 
linked to the modulations of the X-ray flux by blobs of matter or- 
biting in the accretion disk. The upper peak is assumed to be the 
modulation at the azimuthal frequency whereas the lower one 
that at the periastron precession frequency v p = — v$ of the orbit 
of the orbiting blob. 

An observed feature of twin pea k QPOs is the 3:2 ratio at which 
their central frequencies cluster. lAbramowicz & Kluzniakl J200 lh 
proposed that the twin peaks could be excited by resonance mecha- 
nisms between oscillation modes in a curved space-time. Later, the 
idea of the 3:2 f requency cluster cau sed by resonant mechanisms 
was reviewed by Bell oni et alj d2005h . who claimed that a random- 
walk could reproduce the cluster at the 3 :2 ratio. Subsequen tly the 
random-walk hypothesis was debated bv lTorok et alj J2008). 
The shape of the modulations by a hot-spot orbiting in the Kerr 
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1 r g = GM/c 2 is the gravitational radius of the compact object 



2 v r is the oscillation frequency of the radius of the orbit: from the perias- 
tron to apoastron and back. 
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metric was studied in detail by Schni ttman & Bertschingej fc004) 
by means of ray-tracing techniques. The power spectra of the sig- 
nal are characterized by peaks at the azimuthal frequency v<j, the 
beats v^, + v T and their Fourier harmonics. It would be interesting 
to investigate whether the beats + v r originate from the timing of 
the azimuth phase <p(t) of the orbiting object rather than from the 
modulations on the amplitude of the signal because of relativistic 
effects (Doppler/gravitational shifts). Hence we concentrate on a 
frequency m odulated (FM) rather th an on an amplitude modulated 
(AM) signal jTaub & Shillindll986l) . This is one of the issues that 
has stimulated the investigat i on pre sente d in this lette r. 

Works by ICadez et all (2008) and iKostic etafl j2009h high- 
lighted the role of the strong tidal force on matter orbiting close to 
a black hole. These works modeled the N IR/X-ray flares th at are 
observed coming from the galactic center dGenzel et al.l2003|). The 
nume rical code on tidal disruption of small satellites ItCostic et all 
2009) was used to fit the observed light curves. 
The simulations of the tidal stretching of a clump of matter were 
applied to the twin peak HF QPOs as well. It was possible to repro- 
duce the high-frequency part of the typical power spectra observed 
in LMXBs. The simul ated power spectra s how the power law with 
the twin peak QPOs dGermana e t al. 2009). 

It is not yet clear whether the upper/lower peak of the twin peaks 
corresponds to the azimuthal frequency. The results from the tidal 
interaction code show that the lower peak in frequency of the twin 
peaks corresponds to the azimuthal frequency y^, whereas the up- 
per one corresponds to the m odulation at v A + v r . This is different 
from the RPM assumptions ^Stella et alii 19991) . since it links the 
lower peak to v p = - v r and the upper one to v^. 

In this letter we highlight the temporal behavior of the az- 
imuth phase 4>{i) of an orbiting object and suggest it as root of the 
observed twin peak QPOs. Our goal is to investigate the appear- 
ance/disappearance of twin peaks in the power spectrum, inves- 
tigating as well what modulation corresponds to the upper/lower 
peak: v 4 ,,v 4 , + v r or v - v r ? 



2 ORBITAL MOTION ON NON-CLOSED ORBITS 

The repeating pattern in time of the azimuth phase <f>(i) for orbital 
motion on a circular orbit can be described by a simple sinusoid 

A sin (2ntv + <p ) (1) 

where A is the amplitude, (f>o the phase at t = and v = 1/P the 
oscillation frequency of the sinusoid of orbital period P. The linear 
term 2ntv draws the linear increase in time of the azimuth phase 
4>(t), i.e. its timing law is t/>(t) = 2ntv + <po radians. 
On an eccentric orbit the orbital motion is faster at the periastron 
than at apoastron. Thus, the phase of the azimuth on an eccentric 
orbit is not a linear function of time. We shall write the orbital mo- 
tion on an eccentric orbit 

A sin (^.Titv^ + B \r p , ej sin (2ntv r ) + 0o) (2) 

where B(r p , e) is the amplitude of the oscillating phase because 
of the different orbital speed of the object at periastron/apoastron 
passage, happening at the radial frequency v r . We can write 
B(r p , e) ~ 2n(\ p -\ a )/w radians, where v p and v a are the orbital 
speeds at the periastron and apoastron, v is the mean orbital speed 
on the orbit of periastron r p and eccentricity e. The exact value of 
B(r p , e) is calculated numerically in Section |2~Tl 

In a curved space-time the orbits are non-closed, v T is always 
smaller than v^. This leads to the periastron precession of the orbit 
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Figure 1. Top: Phase of the azimuth <p(t) for orbital motion on a non-closed 
orbit in the Schwarzschild metric. Middle: Phase-residuals after subtracting 
to (f>(t) the best-fit first-order polynomial. Bottom: Phase-residuals for an 
orbit at a different radius (see text). 



at the frequency v p = - v r , like in the Mercury's orbit. The power 
spectrum of equation (O might give peaks besides that at v^. 

The observed frequency of the twin peak QPOs are typical of 
the orbital motion in a region of the space- time within r = 10 r g 
from the compact object. For a Schwarzschild black hole the ratio 
v r /v<j ranges from at r = r ISC o = 6 r g (at ISCO v r = 0) to 1 
at infinite r. Thus, from equation (O the phase of the azimuth is a 
linear function of time only at ISCO. 

We write the FM signal 

S (t) = Random (?) + A sin {2ntv,p + B(r p , ej sin (2ntv r ) + O ) (3) 

where the term Random(f) is random noise, as expected in the sig- 
nal from LMXBs. The power spectrum of 5 (r) is the squared mod- 
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ulus of the Fourier transform S'(v) calculated numerically, cross- 
checked with the task powspec in the timing analysis software 
package Xrono j^. 

We should bear in mind that the exact shape 5(f) of the 
modulation by an orbiting object/hot-spot is obtained by per- 
forming ray-tracing of photons in the curved sp ace-time (e.g. 
ISchnittman & Bertschingeill2004l : [Kostic et al]|2009l) . Here we are 
only attempting to study the azimuth phase temporal behavior 
(f>(t) of the orbiting object/hot-spot, i.e. its kinematics. This is a 
well known method used in pulsar timing analysis to study the 
neutron star and its environment. Pulsars timing is performed by 
studying the phase temporal behavior <ft(f) of the rotating hot- 
spot (e.g.lHulse & TavloJl975l:IWolszczajl994lstairs et alhool 
I Kramer et a l. 2006; Germa na et all 20 121) . 



2.1 Orbital motion in the Schwarzschild metric 

To justify the approach described above, we calculate the orbital 
m otion in the Schwar zschild metric, following the parametrization 
in lCutler et alj JT994 Fl. Fig. [Tj (top) shows the parametric plot be- 
tween the azimuth phase <p(x) and the coordinate time t(x), for an 
orbit of eccentricity e = 0.1, with periastron at r p = 7.5 r g and 
around a 8 M Q black hole. The overall linear increase simply shows 
the amount of the azimuth phase in coordinate time, as function of 
the radial phase parameter^. The simulation lasts for 10 radial os- 
cillations, corresponding to 20 azimuth turns (~ 125 rad), because 



for this orbit = 2v r . 

The overall linear behavior in Fig [T] (top) is simply the argu- 
ment 2ntv<p in equation 10. For a circular orbit, only the term Intv^, 
describes the orbital phase 4>{i). We notice that in Fig [T] (top) the 
overall linear behavior wobbles. Fig [T] (middle) shows the residu- 
als of the azimuth phase after fitting to the function <f>{t) a first-order 
polynomial and subtracting it. There are oscillations at the relativis- 
tic radial frequency typical for this orbit, v T = 90 Hz. The amplitude 
of the oscillation is B(7.5 r g ,0. 1) ~ 0.3 rad. Fig [T] (bottom) shows 
the phase-residuals for an orbit with the same e at an inner radius, 
r p = 6 r Jl- At this radius, v r = 70 Hz- In FigQ~]bottom the number 
of oscillations is smaller than in Fig[T]middle, because v, is smaller. 

To conclude, the azimuth phase t/>(t) on a slightly eccentric or- 
bit close to a compact object (r ~ r g ) is not a linear function of 
time. The phase oscillates at the radial frequency v r because of the 
different orbital speed at periastron/apoastron passage. In a curved 
space-time the frequency at which an orbiting object is at the peri- 
astron/apoastron is v r t v$. We write the timing law 



</> (f) = 2ntv^ +B(r p , e) sin (iMv,-) + < 



(4) 



Because the first derivative of the phase <p{i)l(2n) is the frequency 
of the orbital motion as a whole, in a curved space-time the orbital 
motion is not periodic. 



3 http://xronos.gsfc.nasa.gov/ 

The integrals reported in ICutler et alHl994l) . describing both the azimuth 
phase <p(x) and time t(x) as function of the radial phase parameter x, were 
calculated in the limit of orbits of small eccentricity. The relativistic fre- 
quencies from this calculation were cross-checked with those obtained an- 
alytically (Germana 2006, Laurea thesis). 

5 For an orbit with e = 0. 1 the innermost stable bound orbit has the perias- 
tron at r„ = 5.5 r„. 




Figure 2. Top: The simulated FM signal S(t) from a non-closed orbit de- 
scribed by equation (3j- The bin time of the temporal series is dt = 0.0001 s. 
Bottom: Power spectrum of the signal described by equation (5) (B = 0.5). 



3 RESULTS 

The noisy function S(t) (eq. ]3]) is shown in Fig. [2] (top). From 
the calculations in Section [2~T1 the radial frequency is v r = 90 Hz 
whereas the azimuthal one is = 2v r = 180 Hz, for an orbit with 
periastron r p = 7.5 r g around a 8 M Q Schwarzschild black hole. We 
are not interested in any parametrization of the amplitude of S (t), 
so we choose to use arbitrary units to describe it. 



3.1 Power spectra 

Fig. [2] (bottom) shows the power spectrum of the FM signal 
5(r) described in Fig. [2] (top). The simulated power spectrum 
reproduces the twin peaks seen in observed power spectra of 
LMXB XTE J1550-564 contain ing a black hole dOrosz et al .120021; 
iRemillard & McClintockl 120061) . In Fig. |2] (bottom) the upper- 
most couple of peaks in frequency corresponds to = 180 Hz 
and + v r = 270 Hz. Because in this region of the space-time 
(r p = 7.5 r g ) v r /v^ = 0.5, + v, and are in a 3:2 ratio. In ob- 
served power spectra the peaks are in a 3:2 ratio. This feature 
has led to the idea t hat twin peaks could be produ ced by res- 
onance mechanisms dAbramowicz & K luzniak 2001). The simu- 
lation shows a third peak at the periastron precession frequency 
v p = v ip — v r = 90 Hz. The three peaks display a harmonic re- 
lation in frequency, 1: 2:3. Such harmonic relation was noticed by 
IRemillard etail (2002) in a set of observations of the source XTE 
J1550-564, with the frequency of the peaks ~ 90: 180:270 Hz. 

Fig.[3](top) shows the power spectrum of the signal S (f) for an 
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Figure 3. Power spectra of the signal described by equation {3} showing the 
uppermost two peaks in frequency not in a 3:2 ratio. Top: Power spectrum 
with twin peaks for an orbit with r p ~ 10 r g and B = 0.5. Bottom: Power 
spectrum for an orbit at r p ~ 6 r„. 



orbit at r p ~ 10 r g with the uppermost couple of peaks not in a 3:2 
ratio. The azimuthal frequency is = 108 Hz and the radial one 
v r = 76 Hz, such that (v^ + v r )/v^ ~ 1.7. Fig. [3] (bottom) shows the 
power spectrum for an inner orbit, r p ~ 6 r g , with = 238 Hz and 
v r = 70 Hz. The ratio of the uppermost couple is (v^ + v r )/y (4 ~ 1.3. 
Fig.[3](top/bottom) suggests that resonance mechanisms to explain 
the twin peaks might not need to be invoked. 

The ratio (v<j + v,)/v<j ranges from ~ 1.7 to 1.3 as long as 
the frequency ranges from 108 Hz to 238 Hz, thus with a rel- 
ative increment of (238-108)/238~ 0.55. We can use this clue to 
further investigate whether the couple of twin peaks in the ob- 
served power spectra is (v^, + v r ) or (y^ - v r , v^). If the ob- 
served couple of peaks were (v^ - v r , y^) they would have been 
in a 3:2 ratio at r p ~ 6.3 r g . Around such radius (r p ~ 6.5 - 5.8 r g ) 
the ratio ranges fr om 1.7 to 1. 3 as lo ng as relatively increases 
by 0.19. Fig. 3 in lTorok et al.l l|2006) shows the frequency distri- 
bution of twin peaks in several sources. In NS LMXBs the dis- 
tribution clusters at the 3:2 frequency ratio and mostly ranges from 
~ 1 .7 to 1.3 (the lower peak from 400 to 900 Hz, the upper one from 
700 to 1200 Hz). If the observed couple of peaks is (v<j - v,, v^) 
the upper peak on the y-axis of that figure relatively increases 
by 0.42. The typic al relative error on the m easured frequency of 
HF QPOs is ~ 1% ( Ivan der Klis et alJ l 19961). implying an error of 
(cr (V2 _ vl) /(v 2 - vi) + o- ri /v 2 )xOA2 ~ 0.02 on the factor 0.42. Thus, 
the relative increment of 0.42 from observations is statistically in 
disagreement with the numerical one, equal to 0.19. If instead the 



observed couple of peaks co rresponds to {v^ju + v r ) the lower peak 
on the x-axis of Fig. 3 in lTorok et al.l d2006l) relatively increases 
by ~0.56, in agreement with the 0.55 increment estimated numeri- 
cally. This further suggests that the observed twin peaks correspond 
to the couple (v^, + v r ) and they might take place over a window 
of radii corresponding to r p ~ 6 - 10 r g in the Schwarzschild met- 
ric. 



4 DISCUSSION AND CONCLUDING REMARKS 

The investigation reported in this letter shows that the twin peak 
HF QPOs observed in LMXBs might be the clue of the kinematics 
of orbital motion in a strongly curved space-time (r ~ r g ). With 
respect to a static observer at infinity, the azimuth phase (p(t) is 
described by the linear function 2ntv^ plus a term of amplitude 
B(r p ,e) oscillating at the radial frequency v T t v^. The phase 0(f) 
wobbles at the radial frequency v r making the whole orbital mo- 
tion not being periodic. The observed width of the QPOs, therefore 
their quasi-periodicity, could be a consequence of the orbital phase 
wobbling. 

The appearance/disappearance of the twin peaks depends on 
B(r p , e). The function B(r p , e) characterizes the non-linear part of 
the azimuth phase 0(f). The simulated power spectra show that the 
uppermost couple of peaks in frequency corresponds to + v r and 
v^. Three peaks are reproduced by these simulations and the low- 
est peak in frequency corresponds to the periastron precession fre- 
quency of the orbit v p = v$ - v r . At a given radius (r p ~ 7.5 r g ) 
such three peaks display a 1:2:3 harmonic relation in frequency. 
E vidence of these peaks in a 1:2:3 harmonic relation was reported 
in lRemillard et alj | |2002|) . 

Simulations of the tidal disruption of a cl ump of matter by 
a Schwarzschild black hole dKostic et alj 12009) show that a sig- 
nal with on ly the uppermost two peaks can be obtained, i.e. 
and + v r dGermana et alj|2009l) . The tidal stretching might blur 
the phase temporal behavior of the signal, thus affecting the ap- 
pearance of the beats ± v r . The typical power law seen in the 
observed power spectra might originate from the increase of the 
signal because of the tidal disruption of the clump of matter (e.g . 
iBelloni & HasingeJl990l : lGermana et al.l2009l : lKostic et alj2012l) . 
Thus, tidal effects on matter orbiting a compact object might be 
complemen tary to the featu res of FM signals described here. In- 
terestingly, [RejsTetal] d2012l) recently reported a QPO in the X-ray 
flux coming from the tidal disruption of a star by a supermassive 
black hole. 

We note that twin peaks might be produced even when their 
frequencies are not in a perfect 3:2 ratio. This suggests that 
their nature could n ot be a key feature of resonance mechanisms 
(Bell oni et al .120051) . but rather their root might be inherent to a FM 
signal characterized by the orbital phase (f>(t) on a generic eccentric 
relativistic orbit. 
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